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COUNTING COLORED PLANAR MAPS 
FREE-PROBABILISTICALLY 

ABDELMALEK ABDESSELAM AND GREG W. ANDERSON 



Abstract. Our main result is an explicit operator-theoretic formula for the 
number of colored planar maps with a fixed set of stars each of which has a 

^ fixed set of half-edges with fixed coloration. The formula bounds the number 

■^l,^ ' of such colored planar maps well enough to prove convergence near the origin 

of generating functions arising naturally in the matrix model context. Such 
convergence is known but the proof of convergence proceeding by way of our 
main result is relatively simple. Besides Voiculescu's generalization of Wigner's 
semicircle law, our main technical tool is an integration identity representing 
the joint cumulant of several functions of a Gaussian random vector. The latter 
identity in the case of cumulants of order 2 reduces to one well-known as a 
means to prove the Poincare inequality. We derive the identity by combining 

r^ ' the heat equation with the so-called BKAR formula from constructive quantum 

"^^ field theory and rigorous statistical mechanics. 



CN , 1. Introduction 

> 

l/^ ■ Let Perm(n) denote the group of bijective maps of (n) = {1, . . . ,n} to itself. 

OO I Let c{9) denote the number of cycles mto which 9 e Pcrm(n) decomposes. Let 

Match(n) C Perm(n) denote the subset consisting of fixed-point-free elements of 
order 2. Let Color(n) denote the set consisting of maps of (n) to itself. Given any 

CO . 9 G Perm(n) and 7 e Color(n), we define finite sets 



(1) Map(6',7) = <( t e Match(n) 



9 and t generate a group 

acting transitively on (n) ^ and 

and furthermore 704 = 7 

(2) Mapo(6l,7) = {l eMa.p{9,j) \ c{9) + c{9i) ^ 2 + n/2}. 

Our main result, namely Theorem 12. 3. 51 below, gives a novel free-probabilistic rep- 
resentation of the cardinality |MapQ(6',7)|. We will locate our main result more 
precisely with respect to the literature after formulating it in the next section. 

The significance of the quantity |MapQ(0,7)| in random matrix theory is ex- 
plained by the following well-known formula which we owe to the physicists. We 
refer to f34j for background on matrix integrals and maps. See also fT9] for a good 
explanation to mathematicians of the physicists' viewpoint on matrix models. Re- 
call that an N-hy-N GUE matrix S is a random hermitian matrix whose entries 
have a centered Gaussian joint distribution satisfying ES(j, j)S(«', j') — dijidi'j. 
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2 ABDELMALEK ABDESSELAM AND GREG W. ANDERSON 

Theorem 1.1 ('t Hooft ^). Fix 9 £ Pcrni(7i) and 7 e Color(n). Let 

= (ii,i • ■ ■ ii,ni) • ■ ■ [ik,i ■ ■ ■ ik.uk) ("1 H h J^fc = "•) 

be the decomposition of 6 into cycles. Let Sj^ , . . . , Sn be independent N-by-N 
GUE matrices. Then we have 

(3) |Mapo(0,7)l - Jm^ ^ ^2+n/2-fe ^^' 

where k(-) is i/ie jomt cumulant functional. 

Recall that for C- valued random variables Xi, . . . , Xk with absolute moments of all 
orders the joint cumulant k{Xi, . . . , Xk) € C is defined by the relation 



(4) log J2 E(Xr...X,^^)^ 



''1 ' ■ ■ ^fc 



1^1,.. 


■,vk= 


=0 


00 

>; 


K 


( 


■■,'^k-- 


= 1 


\ 



i/i!-- -i^fc! 




Xi , . . . , Xi , 



standing between formal power series. For k(-) one also has an expression (see 
equation ([M)) below) involving the Mobius function of the lattice of set partitions. 

Of course Theorem II. II is only a faint reflection of the full matrix model picture. 
Formula ([3]) merely peels off the leading term of the 1/A^-expansion. We will not 
be considering the higher order terms in this paper — but see Remark 14.2.41 below. 

Physicists have a visually appealing interpretation of t g MapQ(0,7) which we 
recall to justify the title of the paper. One draws a diagram as follows. First one 
marks down some vertices in the plane indexed by the 0-cycles. Then out of each 
vertex one draws half-edges indexed by the elements of the corresponding 0-cycle, 
arranging them in the circular order dictated by 9 and "coloring" them by 7, thus 
forming a "star." Finally, one joins each half-edge to its same-colored mate via the 
perfect matching t to form a whole edge, taking care that connecting paths do not 
cross — the numerical condition c{9) + c{9i) =2 + n/2 guarantees the possibility of 
drawing the picture in the plane without crossings. The resulting graph embedded 
in the plane (up to some abuse of language) is a colored planar map. See (35] for 
background on graphs embedded in surfaces. 

Actually our work in this paper consists entirely of an analysis of the right side of 
equation (jH]). Thus the reader could without loss of comprehension of our arguments 
take ([3]) as the definition of the quantity |Mapo(0,7)|. 

Although we must defer the statement of Theorem 12.3.51 we can immediately 
state a couple of corollaries to it in order to convey its flavor. We use the notation 
Il{- • • } for the characteristic function of the condition between braces. For constant 
7 G Color(n) (the "monochrome" case) we just write Mapo(0) instead of Mapo(0, 7). 

Corollary 1.2. Let ni, . . . , n^ > be integers. Put n ~ ^Ui and p — Y[^i- Let 

9 G Pcrm(n) have cycles of length rii, . . . , n/j. We have 

(5) |Mapo(6')| < pn'''-22«-2fc+2jj|^ > 2fc - 2}. 

This estimate may be new. In any case the method of proof is surely new. See ^2.3.61 
below for proof of the corollary, and see Proposition l2.3.7] immediatelv following for 
a more precise version of this bound. 
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Corollary 1.3. Let ni,...,nk > be integers. For integers vi,...,Vk > let 
Oui,...,Uk G Perm(^ z^irii) have cycles of length ni, ... ,ni, ... ,nr, ... ,nk. Then the 

generating junction 

CO OO Vi Vk 

(6) E---E|MaPo(^.......JI^^r^ 

converges for (zi, . . . , z^) in a neighborhood of the origin in C'^. 

This result is certainly not new. It follows e.g. from results of [T7] or [I^]. But the 
relatively elementary character of our proof is novel. 



Proof. By Corollary [TT2] the series in question is majorized by the series 

and the latter clearly has a positive radius of convergence. D 

Now we turn to a topic which seems at first glance only mildly relevant. Let 
Ci C^"'^'*! C*'^^ G M" be independent random vectors, each with i.i.d. standard Gaussian 
entries. Let /, g : M" — > C be adequately nice functions. The identity 

(7) Cov(/(C),. 9(C)) = I E (V/(C(^)) • Vg (tC,^^^ + yi^^C^'^)) di 

is well-known. Notably, ([7]) implies the Poincare inequality with the best constant. 
For background, further applications of (O and more references, see e.g. |9]. The 
main technical result of the paper, namely Theoreni l2.4.5l below. is a generalization 
of ([7]) holding for joint cumulants of arbitrary order. We will locate our main tech- 
nical result more precisely with respect to the literature after formulating it in the 
next section. We will derive Theorem l2.4.5l from the Brydges-Kennedy-Abdesselam- 
Rivasseau (BKAR) formula [Ml [5] which has appeared in the context of constructive 
quantum field theory and rigorous statistical mechanics. 

Here is the plan of the paper. In f}2]we formulate our main result, Theorem l2.3.5l 
We also state our main technical result. Theorem 12.4.51 In fJ31 we review the 
BKAR formula (see Theorem 13.1.41 below) , give pointers to the related mathe- 
matical physics literature, and finally we use BKAR and the heat equation to 
prove Theorem 12.4.51 In 21 ^^ prove Theorem 12.3.51 by combining Theorem 12.4.51 
Voiculescu's multi-matrix generalization of Wigner's semicircle law [321 Thm. 2.2] 
and Theorem ll.il 



2. Formulation and discussion of the main results 

2.1. Co-ultrametrics, forests and random matrices. We introduce notions 
needed to state our main result. Theorem 12.3.51 below, and our main technical 
result. Theorem 12.4.51 below. As in the introduction, we write (fc) = {1, . . . , fc}. 
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2.1.1. Co-ultrametrics. Let A be a k-hj-k real symmetric matrix whose entries 
A{i^j) satisfy the fonowing conditions for all ii, J2,i3 G {k): 

(8) A(zi,*i) = 1, 

(9) A(ii,z2) e [0,1], 

(10) A(ii,i3) > uAn{A{ii,i2),A{i2,h))- 

We call A a fc-by-fc co-ultrametric. We choose this terminology because (fTO|) is the 
reversal of the ultrametric inequality familiar (say) to number theorists. Let coU(fc) 
denote the set of fc-by-fc co-ultrametrics. 

Lemma 2.1.2. Fix A e coU(fc). Let ji, . . . ,jr G (fc) be a sequence with r > 2. 
Then the two smallest numbers on the list A{ji,J2), ■ ■ ■ , A(jV,Ji) are equal. 

This is the "co" version of a standard fact about ultrametrics. ("Every triangle is 
isosceles.") 

Proof. On the one hand, after cyclically permuting the indices j'l, . . . , j^, we may 
assume that A{ji,jr) = A{jr,ji) < min(A(ji, J2), • • ■ ,^(>-i, jV))- On the other 
hand, the reverse inequality holds by repeated application of ([TU|. D 

Lemma 2.1.3. Every A G coU(fc) is positive semidefinite, and moreover positive 
definite if all off-diagonal entries of A are less than 1. 

Proof. For every A S coU(fc) and t e [0, 1) there exists unique A^^*^ € coU(fc) 
with entries A'^*\i,j) = l{A{i,j) > t}. One checks immediately that A*^*^ is the 
indicator of the graph of an equivalence relation in (fc) and in particular is positive 
semidefinite. Furthermore, if the off-diagonal entries of A are all less than one, then 
A^*) is the identity matrix for t near 1. The trivial formula A = L A^*) dt completes 
the proof. D 

2.1.4. Gapless co-ultrametrics. We say that A G coU(fc) has articulation equal to 
the cardinality of the set {A{i, j) \ i,j = 1, . . . ,k}\ {0, 1}. We say that A £ coU(fc) 
has CO- articulation equal to the number of equivalence classes in the set (fc) for the 
equivalence relation i ^ j <i^ A{i,j) > 0. We say that A S coU(fc) is gapless if 
articulation and co-articulation sum to k. 

2.1.5. Graphs (especially trees). In this paper a graph — (V, E) is a pair consisting 
of (i) a finite set V of vertices and (ii) a set E of (unoriented) edges each of which 
is a two-element subset of V. We say that © = {V, E) is a tree if © is connected 
and \E\ = \V\ - 1. We say that a tree 1 = {V, E) spans a set 5 if 5 = V. RecaU the 
generating function identity 

(11) Y. n ^'^^- = (^i+-' ' -'"-^ 

'X: tree {'^^j}'- edge of T 

spanning {k) 



Xk) Xi---Xk 



which one obtains by specializing KirchhofF's matrix-tree theorem. In particular, 
there are exactly fc*^"^ trees spanning the set (fc). The graphs of interest in this 
paper are mostly trees and if not trees then forests, i.e., disjoint unions of trees, 
or equivalently, graphs with no circuits. When we say that 5' is a forest in (fc) we 
mean that ^ = {V,E) \s a. forest with vertex set V = (fc). 

In fact the concept of gaplessness is closely related to forests, as follows. 
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Proposition 2.1.6. Let ^ = {V,E) be forest in (k). Let {xe}eeE be a family of 
distinct numbers selected from the open unit interval (0, 1). Then there exists unique 
A e coU(A:) such that 

(12) A{i,j) = X[ijj for every edge {i,j} G E and 

(13) A is gapless of articulation \E\. 

Proof. By Lemma 12.1.21 for distinct i^j e {k) connected by some walk in 5", any 
co-ultrametric A satisfying (fT2|) and ()13|) also satisfies 

(14) A(i,j) =min(a;{j-^j^},...,a;{j-_^j„}) 

where j'l • • • jr is the unique geodesic walk in "^ from i to j. (It is at this point that 
the assumption of distinctness of the numbers x,, enters crucially.) Now consider 
the following two equivalence relations ~i and ^^2 in (k). Let i ~i j if and only if 
i and j are joined by a walk in ^. Let i ~2 J if and only if A{i,j) > 0. We have 
* ~i J =^ * ~2 j- But furthermore, ~i has k — \E\ equivalence classes because ^ 
is a forest and ~2 has k — \E\ equivalence classes by (|13p . Thus the equivalence 
relations ~i and ~2 coincide. It follows that 

(15) A{i,j) = if i, j e (fc) are not joined by a walk in 5^. 

Thus uniqueness is settled. Conditions ([5]), (|14p and (IT5|) define a fc-by-fc real sym- 
metric matrix which one can easily verify does indeed belong to coU(fc) and satisfy 
p2|) and (|T3)) . Thus existence is settled. D 

The preceding result has an easily derived converse. 

Proposition 2.1.7. Let A £ coU(fc) be gapless of articulation £. Let 

{iiJi},---,{iiJi} 
be two-element subsets of (k) such that 

0<A{iiJi) <■■■ <A{ii,,je) < 1. 

Then the graph ((fc), {{ii,ji}, • ■ • , {ie,je}}) is a forest in (k) . 

Proof. Lemma [2 . 1 . 2 1 rules out the possibility of circuits. D 

2.1.8. The weight matrix attached to a forest. Let ^ = {V,E) be a forest in (fc). 
Let {Xe}eeE be i.i.d. random variables uniformly distributed in (0,1). We define 
wtj G coU(fc) to be the random matrix which with probability 1 is uniquely deter- 
mined by Proposition 12 . 1 .61 By Lemma 12.1.31 with probability 1, the matrix wty 
is positive definite and thus has a unique positive definite square root -^wtj. We 
introduce the random matrices wty just to have a convenient compact notation 
with which to handle integrals on various pieces of the space coU(fc). 

2.2. Spanning trees and splicing involutions. We introduce some partly group- 
theoretical and partly arboreal notions needed to formulate Theorem 12.3.51 below. 

2.2.1. Permutations. Recall from the introduction that Perm(n) denotes the group 
of bijective maps of (n) = {1, . . . , n} to itself and c{9) denotes the number of cycles 
into which 9 e Perm(n) decomposes. We call r e Perni(n) an involution if t^ = 1. 

Proposition 2.2.2 (Riemann-Hurwitz bound). Let ao^ai,aoo € Perm(n) generate 
a group of permutations acting transitively on (n) and satisfy O'oo'ifoo = 1- Then 
we have c{ao) + c(cri) + 0(0-00) < n + 2. 
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Proof. The triple (co, ctijCToo) is an example of a constellation. See [22j for back- 
ground on constellations. To the triple (cto , fi , CToo ) one naturally attaches an 
n-sheeted covering of the Riemann sphere branched at 0, 1 and (X). The genus 
g of that covering, i.e., the number of its "handles," satisfies the Riemann- Hurwitz 
formula 2g — 2 — ~2n + J2i=o i oo('^~ c((Ti)). The result follows because g >0. D 

Remark 2.2.3. Note that the numerical condition figuring in the definition ^ of 
Mapo(0,7) touches the bound enunciated in Proposition 12.2.21 and thus implies 
"planarity" of elements of Mapg(6', 7). 

2.2.4. Splicing involutions. Let v : (n) — > (fc) be an onto function. Let T be a tree 
spanning (fc). We define the set of splicing involutions 

Splicefj(j/) C Perm(n) 

indexed by T and v to be the subset consisting of involutions t admitting a factor- 
ization 

T= («i,ji)---(«fe-iife-i) 

into disjoint transpositions such that 

1 = mAMn), Hn)}, • ■ • , M^k-i), Hjk-i)}})- 

Proposition 2.2.5. Let 9 G Perm(n) be any permutation and put k — c{9). Let 
V : (n) — > (fc) be any 0-invariant onto map. Let T be any tree spanning (k). Then 
for all T € Splicefj(z^) the composite permutation 6t is cyclic. 

Proof. Since onto, the map ly in effect enumerates the 0-cycles. It follows in turn 
that the disjoint transpositions into which t factors provide enough linkages between 
the 6'-cycles to make the Cayley graph of the pair {9, t) connected. Thus 9 and r 
generate a group of permutations acting transitively on (n) . By Proposition 12.2.21 
it follows that c{9t) = 1. D 

Proposition 2.2.6. Let 1/ : (n) —¥ (k) be onto. We have 

na^ \- iq r ^^1 fi^i^ffili^ ifn>2k^2, 

(16) 2^ |Splice2:(z/)| = <^ {n- 2k + 2)1 '' " 

1: tree I ifn<2k-2, 

spanning (k) 

where ni = \h'^^{i)\ for i = 1, . . . , fc. 

Proof. For each tree T spanning (fc), let T be the set of ordered pairs 
{i,i') € (fc)^ such that {«,«'} is an edge of 1, let tt^ = ((i,i') i->- i) : T — )■ (fc), 
and put fii^i = JTT^ (i)| for i = 1, . . . , fc. Now the set Splice<j(z^) is in evident bi- 
jective correspondence with the set of one-to-one maps ip : 1 -^ (n) such that 
i> o ip = 771. Thus we have 

k 

{LRSoim)= E U 



%: tree 
spanning (fc) 






II' 



where it is understood that we set , /^''^ -,, equal to zero if fi-^^i > Hi. Note 
furthermore that we can rewrite identity (jlip in the form 

k 

X] n ^'^•i(<=') = X! n^i'^' ^ (^1 ^ — hxk)''~^xi---xk. 

T: tree e^X ^' ti"cc «— 1 

spanning (k) spanning (k) 
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Finally, we have 

which finishes the proof. D 

2.2.7. Colored splicing involutions. Let i^ : (n) — > (A;) be an onto function. Let T 
be a tree spanning (fc). Given also 7 e Color(n) we put 

Splice<j(i', 7) = {r G Splice2;(z^) | 7 o r = 7}, 

thus defining the colored version of Splice<r(i^). 

2.3. Statement of the main result. 

2.3.1. Data. Fix 6 S Perni(?T,) and 7 E Color(n). Our aim is to formulate a result 
representing the quantity |Mapo(^,7)| free-probabihstically. For this purpose we 
put k = c{6) to abbreviate and we fix a ^-invariant onto map v : (n) — > (fc). 

2.3.2. The canonical splicing polynomial. Let C({Xi}"^]^) be the noncommutative 
polynomial algebra generated by a family {Xi}"^]^ of n independent noncommuting 
variables. Let T be any tree spanning (fc). We then define 

rGSplice.j-(i/,7) ii ,...,in€{n) 

s.t. 9T—{ii---in) 
and ii — 1 

The inner sum has exactly one term by Proposition 12.2.51 so that the number of 
monomials in Polyg ,j,^ .j equals |Splice3;(;^, 7)|. Clearly, Poly^,^ ,^<j: is homogeneous 
of degree n — 2fc + 2 in Xi, . . . ,X„ and in particular vanishes identically unless 
n> 2k — 2. One should think of Poly^ ^r ^^ a sort of noncommutative generating 
function representing the set Splice5;(z^, 7). 

2.3.3. A transformation of the canonical splicing polynomial. Let {{'Zi{i, j)}^^i}'^^i 
be another family of independent noncommutative algebraic variables. Given a tree 
T spanning (fc) put 

k 

Z5(i,j) = ^V^(^0Z(«',j) for(z,j)e(fc)x(n), 
(which is random) and in turn put 



Aihg^^^i, 



E Efpoly,_^,,_j 

■CCS T \ 



trees 
spanning (k) 



for i ^ (n) / 



where the expectation is computed term- by-term, i.e., one first expands the inte- 
grand as a sum of finitely many monomials in the variables 7i{i,j) with bounded 
random coefficients, and then takes the expectation of each of the finitely many 
nonzero coefficients. We remark that Arbe^-,,^^ actually is independent of v. This 
fact will not be needed in the sequel. 



8 ABDELMALEK ABDESSELAM AND GREG W. ANDERSON 

2.3.4. Free standard semicircular variables. Let (.4, Lp) be a faithful tracial 
C* -probability space. Let {{z(i, j)}*L]^}"^]^ be a family of free standard semicir- 
cular variables in [A, (p). A quick introduction to free probability theory more than 
adequate for our purposes here can be found in ^ Chap. 5]. For discussion at length 
and in depth see |3S] or [55] , 

We are ready to state the main result of the paper. 

Theorem 2.3.5. Notation and assumptions are as above. We have 

(17) |Mapo(6',7)| =(^ Arbe,^,^ z(».j)=z(jj) 

V for{i,])e{k)x{n)/ 

The proof of Theorem 12.3.51 will be completed in 21 below. 



2.3.6. Proof of Corollary \1.2[ The corollary is entirely subsumed by the following 
more precise technical result. 

Proposition 2.3.7. In the setup for Theorem \2.3.5\ and with \\-\\ denoting the norm 
on the C* -algebra A, we have 



(18) 



Arbfl 



,7,1^ 



for (i,j) e (fc)x(n) 



yfi-2fc+2 ' 



fc)!nti^ 



< 



2A: + 2)! 



ifn>2k-2, 
ifn<2k~2. 



Proof. Given any tree T spanning (fc) we define 
fc 
Zi{i,j) = 2_^ ■\/wtJ(i,i') z{i ,j) £ A for i e (fc) and j G 



(which is random). With probability one, because the diagonal entries of the matrix 
wt-j: are identically equal to 1, for each fixed i G (fc) the family {z'x(i,i)}?=i is free 
standard semicircular and in particular ||z'j(J, j)|| — 2 for all i and j. Consequently 
we have 



Poly^ 



iT^i^/i 



X,=z.j(i/(j),7(j)) 
for i£ (n) 



< 2 



n-2k+2 



|Splice<j;(i^,7)|, 



almost surely. Trivially, for a tree T spanning (fc) and any bounded linear functional 
lp € A* we have the formula 



„k-2 



EV- Polye 



i7)'^/-i- 



for iG{n) . 



= ^ Arb6 



,7,1^ 



Z(ij')=z(ij) 

for (j,j) e (fc)x(n>^ 



By Hahn-Banach we can choose tp of norm 1 to make the right side above equal to 
the left side of (|18p . We then get the result by Proposition 12.2.61 D 

Remark 2.3.8. If fc = 1 and 6 = {ii ■ ■ ■ in), formula (IT71) reduces to the statement 

(19) |Mapo(0,7)l = ¥'(z(l,7(n)) • • • z(l,7(z„))), 

which nowadays, free probability theory taken for granted, one interprets as an 
instance of the semicircular analogue of the Wick formula. 
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Remark 2.3.9. In a continuation of the previous remark, we note that when spe- 
ciahzed to the case fc = 1, Theorems 11.11 and 12.3.51 in tandem are equivalent to 
Voiculescu's resuh P], Thm. 2.2]. 

Remark 2.3.10. We turn now to the case fc = 2. In this case Theorem [533] provides 
a formula for the same quantity as does [24i Thm. 5.3]. The formulas look rather 
different. On the one hand, [24l Thm. 5.3] involves an extension of free probability 
theory to the second order (cyclic Fock spaces and annular pair partitions inter- 
vene). On the other hand, Theorem l2.3.5l is phrased in terms of the usual first order 
theory, albeit with certain operations of differentiation and integration added on. 
Detailed relations between the two types of formulas remain to be worked out. 

Remark 2.3.11. In the case fc > 2 there is some overlap, at least thematically, be- 
tween our work here and that in [16]. We note, for example, that both theories 
have in common a key role for permutations subject to what we call here Riemann- 
Hurwitz bounds. But the program of [16^ is much more ambitious, seeking to under- 
stand limiting behavior of joint cumulants of traces of monomials in independent 
unitary invariant random matrices (not just GUE matrices), and it does so by en- 
larging the foundations of free probability to all orders. More precise connections 
to the less sophisticated theory worked out here remain to be elucidated. 

Remark 2.3.12. One can construct noncommutative (a priori possibly signed) mea- 
sures through consideration of generating functions of the form ([6|) for small values 
of the parameters. The formidable analyses of [19 and "2^ taken together show 
(among many other things) that under the evident necessary conditions those mea- 
sures are positive and engender von Neumann algebras with useful and interesting 
properties. An important motivation for our work was to find a relatively elemen- 
tary approach to the positivity phenomenon independent of matrix models. We did 
not succeed in finding it but we nonetheless hope Theorem [233] could provide some 
clues in this direction. 

Remark 2.3.13. Tutte developed over several decades and in many papers an ap- 
proach to the enumeration of planar maps based on "well-labeled trees" and gener- 
ating functions. For just one influential example, see [31]. We also mention [8], [10] 
and \TD as recent papers in Tutte's line of descent. The relationship between the 
trees appearing here and those of Tutte is apparently not very direct but deserves 
investigation. 

2.4. Formulation of the main technical result. We state the main technical 
result of the paper, namely Theorem 12.4.51 below. 

2.4.1. The m,atrix space Matfexn- Let Mat^xn denote the space of fc-by-n matrices 
with real entries. We equip Matfexn with the inner product (a;, y) — x ■ y = tY (x'^y) 
and we put l|a;]| — {x,xy^^. We write Matfe — Matfexfe- Let Ife G Matfe denote 
the fc-by-fc identity matrix. Let Dij denote differentiation of smooth functions on 
Matfexn with respect to the entry in position (?, j). Given a matrix Z e Matfexn, we 
often denote its entry in row i and column j by Z{i,j). We tend to choose notation 
of the latter type when there are further indices to keep track of. 

2.4.2. Classes of functions defined on Matfexn- We say that / : Matfexn — ^ C is 
sm,ooth if infinitely differentiable; polynomial if expressible as a polynomial in the 
entries; of exponential growth if there exist constants ci and C2 depending only on / 
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such that |/(x)| < cie*^^"^" for aU x G Matfexn; and tame if smooth and furthermore 
partial derivatives of all orders have exponential growth. We carry the preceding 
terminology over to functions / : R" — >► C via the identification ]R" = Matixn- 

2.4.3. Tensor products. Let /i,...,//£ : M" — > C be tame functions. The tame 
function / : Matfexn ~^ C defined by the formula 



/ 



xii 



a;ir: 



will be denoted hy Ji® ■ ■ ■ ® fk- 



I I /A-^il 5 ■ ■ • ; '^in } 



2.4.4. Differential operators indexed by trees. For each tree T spanning the set {k) 
and smooth function / : Matfexn — > C we define 



^{i.i'}: edge of 'I j = l 

Theorem 2.4.5. Let C, G M" and Z G Matfexn have i.i.d. standard Gaussian 
entries. Assume that for all trees T spanning (fc) the random matrix wtj is inde- 
pendent of Z . Then for any tame functions /i, . . . , /^ : M" — > C we have 



(20) «:(/i(C),...,/fe(C)) 



E E(L<r(/i 



trees % 
spanning {k) 



where k(-) is the joint cumulant functional. 



The proof of Thcorcm l2.4.5l will be completed in §3.4l below after we have introduced 
the BKAR formula. 



Example 2.4.6. In the case fc = 2 of Theorem 12.4.51 consider the unique tree T 
spanning the set (2). Let C/ be a random variable uniformly distributed in (0, 1) 
and independent of the matrix Z . We then have 



wt<j = 



1 U 

U 1 



and 



1 



Z. 



Thus (l^ni) in the case k — 2 reduces to formula ([7]). 



Remark 2.4.7. We will prove Theorem 12.4.51 by exploiting the so-called BKAR 
formula [131 [5] , thus placing it in the long tradition of "arboreal" representations of 
truncated correlation functions in statistical mechanics and quantum field theory. 
But Theorem 12. 4. 51 is also in many respects similar, say, to 26, Thm. 5.1] and thus 
in principle also part of a long tradition of "arboreal" representations of cumulants 
of Wiener processes. To set up a "dictionary" linking these two distinct traditions 
using the specific example of Theorem l2.4.5l is an interesting problem, if largely one 
of exposition. Perhaps the vantage point of the paper [291 permits a unified view. 



Remark 2.4.8. Theorem 12.4.51 complements the infiuential central limit theorem 
|15| Thm. 4.2] by supplying explicit means to compute limiting covariance as well 
as higher corrections in typical random matrix applications. 
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Remark 2.4.9. On the one hand, Theorem 12.4.51 uses more general test-functions 
than we need to make our apphcation to counting colored planar maps. For the 
latter purpose only polynomial test-functions are actually needed. On the other 
hand, the regularity assumptions in Theorem 12.4.51 are clearly far stronger than 
necessary. The minimal hypotheses concerning regularity remain to be worked out. 

Remark 2.4.10. It would be interesting to obtain a cumulant representation for func- 
tions of a vector uniformly distributed on the discrete cube generalizing 
O Thm. 3.1] in the same sense that Theorcm l2.4.5l generalizes ([7]). 



Remark 2.4.11. It would be interesting to find an analogue of Theorem 12.4.51 with 
the Gaussian random vector replaced by a Haar-distributcd element of the iV-by-iV 
unitary group. This problem is suggested by Remark l2.3.11l above. 

3. Application of the BKAR formula 

In this section we state the BKAR formula, derive a "connected" version of the 
formula, provide some background and references, and finally prove Theorem l2.4.5l 

3.1. Statement of the BKAR formula. The BKAR formula is a combinato- 
rial amplification of the Fundamental Theorem of Calculus. In its present form the 
identity first appeared in W, Thm. III.l], building on the previous work of Brydges 
and Kennedy [13^ . For the reader who is not familiar with the literature of construc- 
tive quantum field theory and rigorous statistical mechanics, the most accessible 
reference for BKAR is the recent introduction with complete proofs given in [T]. 

3.1.1. Note on notation. Notation used here is not exactly the same as in p] but 
the differences are of a trivial nature. For example we employ k-hy-k symmetric 
matrices here to play the role assigned in [I] to arrays indexed by two-element 
subsets of (fc). The reader should not have great difficulty translating. 

3.1.2. Notation related to symmetric m,atrices. Let Syruj. denote the space of 
fc-by-fc matrices with real entries. Let Sym^ C Synij. denote the closed set of pos- 
itive semidefinite matrices and let Pos^ C Synij. denote the open set of positive 
definite matrices. For any smooth function / : L/ — >■ C defined on an open subset 
U C Symj,, and given matrices X Cz U and Y e Synij;,, we define the directional 
derivative by Vy/(X) = fj{X + tY)^^^. 

3.1.3. Differential operators indexed by forests. Let ^ — {V, E) be a forest in (fc). 
For any smooth function /:[/—> C defined on an open subset U C Sym^. we define 

(21) ds/ = I n ve..+«.. / 

where eij S Matfc denotes the k-hy-k elementary matrix with entry 1 in position 
(i, j) and elsewhere. If i<^ = 0, then dj is simply the identity operator. 

Theorem 3.1.4 (BKAR formula). For any smooth function f : Sym^, — > C, 

1 ... 1 

' = Y. E(d5/)(wtg). 



(22) / 



forests ^ 
in (fc) 



1 ... 1 

Proofs of Theorem Em can be found in ^ Thm. III.l], [H Thm. VIII.2] and [J. 
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3.2. The connected BKAR formula. To prove Thoorcm l2.3.5l we need to rewrite 
equation ([22]) so as to make the right side into a sum over trees rather than forests. 
We undertake that exercise in this subsection after recahing a few facts about set 
partitions and introducing appropriate notation. See for example [28) or |30) for 
background on the topic of set partitions, Mobius functions and joint cumulants. 



3.2.1. The lattice Part(fc). Let Part(fc) denote the lattice of set partitions of (fe). By 
definition the elements of Part(fc) are disjoint families of nonempty subsets of (k) 
whose union is (k) . Members of a set partition are called its blocks. Given partitions 
$, ^ S Part(fc) we say that <& is a refinement of '^ and write $ < ^ if and only 
if for all A e ^ there exists B G ^ such that A C B. The set Part(/c) is partially 
ordered by the relation of refinement and thus becomes a lattice. Let 

U. = {{l,...,fc}} and Ofe = {{l},...,{fc}}, 

which are the maximal and minimal elements of Part(fc), respectively, with respect 
to the refinement partial order. 

3.2.2. Matrix representation of partitions. Given $ G Part(A:) we define a matrix 
[$] e coU(A:) by the rule that [$](i, j) = Y,Ae<i> Mi''- J) <S ^ x M- For example 



[{{1},{2,3}}] = 



1 1 
1 1 



Note that in general 



[U] = 



and [Ok] 



It is clear that the map ($ i-> [$]) : Part(A;) — > coU(A;) is one-to-one. More precisely, 
for any partitions $, $ G Part(/c) we have $ < ^ if and only if [$](z, j) < [^](i, j) 
for i,j ^l,...,n. 

Remark 3.2.3. It can be shown that for every A G colJ{k) there exist 

• a unique chain {$0 > • • • > ^e} in Part(A;) and 

• unique numbers < ti < ■ ■ ■ < ti < 1 

such that 



^^[$,] + ^t.(['i>.-i] -[$,])■ 



One proves this by elaborating the idea of the proof of Lemma [2. 1.31 Thus coU(A:) 
has a canonical simplicial decomposition indexed by chains in Part(fc) and in this 
qualified sense the points of coU(A:) of the form [$] for <i> G Part(A:) are its vertices. 
But note that coU(A:) is not convex for fc > 3 nor is it a topological manifold. It 
has in general a rather complicated and interesting shape. In the simplest nontrivial 
case fc = 3 it has the fish-like shape of three triangular fins stuck onto a line segment 
backbone. 
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3.2.4. The Mobius function of the lattice of set partitions. The Mobius function 

^{. : .) = ^^(. : .) : Part(fc) x Part(fc) -> Z 
is by definition characterized by the Mobius inversion formula 

(23) Y^ ]i{n < 6} M(e : E) = ]i{n = s} = Y^ ^(n : e) ]i{e < s} 

eePart(5') eePart(s) 

holding for all 11, E e Part(fc). For C- valued random variables Xi, . . . ,Xk with 
absolute moments of all orders it is well-known that 

(24) K{X^,...,Xk)^ J2 MH: 1^)11^11^- 

nGPart(fc) Aen ieA 

This should be compared with formula @ above. 

Corollary 3.2.5 (The connected BKAR formula). In the setup of Theorem \3.1.4\ 
we have 

(25) Y M(n:lfc)/([n])= Y E(dx/)(wt.x). 

nePart(fc) trees 1 

spanning (k) 

Proof Theorem 13 . 1 .41 has the trivial generalization 

/([n])= Y E(d5/)(wt5) fornePart(fc), 

forests 5 = {V, E) in (fc) 

s.t. every edge {i, j} € E is 

contained in some block of IT 

whence the result by Mobius inversion. D 

Remark 3.2.6. In the setup of Corollarv 13.2.51 it can be shown that 
J2 E(d.x/)(wt.x) 

trees T 
spanning (k) 

E , , 

maximal chains "^ "» 0-io<tl<-<tfc-i <tfe-l 

lfe=*o>->*fe-l=Ofc 
in Part(/c) 

Qk-l / fe-1 



/ [^fc-l] + V t,([$,_i] - [$,]) dti • • • dtfe^i. 



dti ■ ■ ■ dtk-i' \^"" "' ' -^ , 

\ *— 1 / 

The proof is a straightforward if tedious application of Propositions 12 . 1 . 6l and [2T771 
along with Remark 13.2.31 We note further that it is not too hard to verify (|25|) di- 
rectly by analyzing the Fundamental-Theorem-of-Calculus-induced relations stand- 
ing between the integrals of the form appearing on the right side above and more 
general integrals of the same form attached to nonmaximal chains in Part(fc). 

Remark 3.2.7. By a straightforward "cutoff" argument, one can allow in equation 
P5p any smooth function / defined on an open set containing Sym^. We will use 
this freedom when proving Theorem 12.4.51 

3.3. Background on BKAR. We provide pointers to the mathematical physics 
literature. 
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3.3.1. Sample applications in physics. While quite elementary, Theorem 13 . 1 .41 can 
be a very powerful tool when aptly deployed. It already has a wide range of appli- 
cations. Here is a non-exhaustive list of such: 

• Explicit convergent expansions for the ground state energy of the spin- 
Boson model [5] as well as Nelson and polaron type models |3] . 

• Uniform L^-type estimates for cumulants (truncated correlation functions) 
in continuous spin models [3]. 

• A close analog of the Ercolani-McLaughlin Theorem in the case of quartic 
interactions p7] . 

• The construction of the Levy area for fractional Brownian motion with low 
Hurst exponent |23j . 

In almost every application, the only difhculty is to "see" the quantity under study 
as / ([Ifc]) for a suitable function /. 

3.3.2. Cluster expansions. The typical application of Theorem 13.1.41 is to provide 
an "arboreal" expression for the joint cumulant of random variables Xi , . . . , X]^ 
with a joint probability distribution (say) of the form exp(— X]q=i ^{xa))AlJ'c{x) 
where V is some function called a potential and d/zc(a;) is a Gaussian measure with 
covariance matrix C. Coefficients s{i,j) are then introduced as multipliers of the 
matrix elements C{i,j) and can be viewed as the "coupling" between the random 
variables Xi and Xj. The game is then to interpolate between the fully coupled 
situation where s{i,j) = 1 and the fully decoupled situation where s{i,j) = 5ij, 
i.e., the situation where the random variables become independent. The resulting 
decoupling expansion generates partially decoupled intermediate terms correspond- 
ing to partitions H of (fc) where random variables belonging to different blocks of H 
become independent. This procedure is called a cluster expansion in the construc- 
tive quantum field theory literature. The first such expansion appeared in |18] . A 
simpler expansion was later introduced in }12j . Theorem 13.1.41 belongs to a third 
generation of yet simpler cluster expansion formulas. 

3.3.3. Generalizations of BKAR. The utility of the representation of cumulants 
provided by the BKAR formula can be assessed by the quality of the bounds for 
such cumulants that one can deduce from it. An important feature is that such 
expansions typically are in terms of trees instead of the much more numerous 
connected graphs with arbitrarily many circuits. Indeed, BKAR produces a sum 
over minimally connected graphs. Actually Theorem 13.1.41 is only one of the sim- 
plest examples of a much more general hierarchy of combinatorial identities which 
can generate hypergraphs instead of graphs and/or graphs which are (say) mini- 
mally p-edge-connected. (The latter corresponds in quantum field theory jargon to 
(p — l)-particle irreducibility.) This very general framework was developed in [6j. 

3.4. Proof of Theorem 12.4.51 We begin with a couple of lemmas, both of which 
are more or less standard. 

3.4.1. The functional %. Given a tame function / : Mat^xn ^ C, we define a 
function Tif : Sym^ — > C by the formula 'Hf{Q) — 'Eif{y/QZ) where as in the 
statement of Theorem 12.4.51 the matrix Z g Mat^xn is random with i.i.d. standard 
Gaussian entries. Dominated convergence and continuity of the matrix square-root 
function implies that Hf is continuous. 
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Lemma 3.4.2 (The independent copies trick). In the setup for Theorem \2.4.5\ we 
have KihiC), ..., fkiO) = E*ePart(fc) ^^i'^ ■ lk){Hf)m) ^here f ^ f, (g> ■ ■ ■ (g> f^. 

Proof. Let {{Cj }?=i}07tAc(fc> be an i.i.d. family of standard Gaussian random vari- 
ables. For $ g Part(fc) and i S (fc), let $(«)€$ be the block to which i belongs. For 
<i> G Part(fc) let Z* £ Mat^xn be the random matrix with entries Z'^{i,j) = (■ . 
For each $ G Part(fc) the law of [$]^/^Z is the same as that of Z*. A straightfor- 
ward calculation then shows that E(/i » • • • (g) /fe)([$]^/^Z) = Has* ^IlieA MO^ 
whence the result by (p4)) and the definition of "H . D 

3.4.3. The functional T-L. Let T-Lf = Hf \posk- Dominated convergence shows that 
Hf : PoSfc — > C is once continuously differentiable, and the latter observation can 
be considerably improved as follows. 

Lemma 3.4.4 (The heat equation). For tame functions f : Mat^xn ~> C and 
i,i' S (k) we have 'Y^^^i'^^ij^i'jf — ^ e^i+eii^f ■ Consequently Hf is infinitely 
differentiable and partial derivatives of Hf of all orders extend continuously from 
Posfc to Sym^ . 

Proof. This is proved by writing 

. . f , ,exp(-itrXTQ-iX)-,^ " 

and taking advantage of the well-known system of partial differential equations 



- Y.^^,0 



o — 1 

acts on entries of Q _ . 

V acts on entries of X/ 



exp(-itrXTg-iX) _ 
(27r)"'=^/2(detQ)"/2 " °' 



We can safely omit further details. D 

3.4.5. Concluding calculation. In the setup for Theorem 12.4.51 where we write 
/ = f\® ■ ■ ■ ® fk and let ^.j- denote summation over all trees T spanning (fc), 
we have 

5]e(Lx/)(V^^) = ^E(HL5/)(wt5) - ^E(d.iH/)(wtT) 
= Imi^ E(d.jH/)(wtT + elfc) = hm ^ ^i{^ : lfc)(H/)([$] + elk) 

% "^ *ePart(fc) 

= Y. M'J':U)(H/)([$]) -^(/i(C),...,/.(C)). 

'I'GPart(fe) 

The first step is justified by Fubini's theorem, the second by Lemma 13.4.41 the 
third by Lemma [3.4.41 and dominated convergence, the fourth by connected BKAR 
fCorollarv l3.2.5p . the fifth by continuity of "H/, and the last by Lemma [3.4.21 D 

4. Proof of the main result 

4.1. The splicing identity. We execute the derivative calculation needed to apply 
Theorem 12.4.51 to the analysis of the expression under the limit sign on the right 
side of ©. 
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4.1.1. Setup. We fix e Perni(n), 7 G Color(n) and an onto 0-invariant map 
v : (n) ~> (fc) as in the setup for Theoreni l2.3.5l We put k ~ c{9). In addition we fix 
a tree T spanning (fc) and a positive integer N. As in the statement of Theorem ll.il 
we fix a decomposition 

S = («i,i • • • «i,ni) • • • («fc,i • • • ik,nj ("1 H \-nk= n) 

into cycles. After relabeling if necessary, we may and we do assume that v{ia,i3) — ot 
for (a, /3) such that 1 < a <k and 1 < /3 < Uq,. 

4.1.2. Commutative variables and a differential operator. Let 

be a family of independent commutative algebraic variables. Consider also the dif- 
ferential operator 



(N) 



n 


n N 


edges {i, i'} 
of 1 


i=l a,^= 



92 



L 



acting on the commutative polynomial algebra generated by the variables Xij (a, j3) 
over the complex numbers. 

4.1.3. Matrices with commutative polynomial entries. Put 



Xij(a,a) if a = /3, 



,(a,/3) = 



Xij(a,fi)+iXij(fi,a) 
Xij{l3,a)-iXij{a,P) 



72 



if a < /3, 
if a>/3 



and let 






Xij{l,l) 



Xy(iV,l) 



x.,(l,iV) 



%(iV,iV) 



for i = 1 , . . . , fc and j = 1 , . 



, ,n. 



We remark that if one evaluates the algebraic variables Xij (a, p) at i.i.d. stan- 
dard Gaussian random variables, then the matrices X^ ■ evaluate to independent 
A^-by-A^ GUE matrices. However, we do not yet carry out that evaluation. We are 
concerned for the moment only with differentiation of polynomial expressions. 

Proposition 4.1.4. Notation and assumptions are as above. We have a relation 



(26) 4"^ {^<^Zl., 

= tr Polye ,^^ 5; 



X 



l,7(ii,«i)^ 



• tra 



(JV) 



■X 



(N) 



)) 



^^l7lyM f°"^ <">/ 



This proposition is the sole motivation for the definition of Poly^ ^ ^ 1. 

Proof. We work in a setup similar to that of the proof of Proposition 12.2.61 Let 

T = {{i,i')e{ky\{i,i'}: edge oil}, 

TT = ((i,j')^j):T^(fc), 
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After opening the brackets in the definition of i^ and noting that 

32 a^ 



dxij {a, I3)dxi'j (a, (3) dxij (a, P)dxi'j (/9, a) 
we get the formula 

AN) ^ V V TT - 

^ ..f M^ ■■! ^ i\ dx^i^^imimMpm- 

k:'i— i-(iV) j:'X-»(n) eel 
s.t.jop=j 

By similarly opening the brackets we obtain identities 

i:{n)^{N} se{n) 



tr Poly9^^,!.^2: 



Xi = X<"^ ,., for i G (ti>. 



_ V^ V^ TT J ^i'M,7M(i(s)>i(^T(s))) if t(s)==s, 

Gsr r ^ wtT/A.. JAl ]l{i(.s) = i(er(.))} if r(.) ^ .. 

TgSphcc,j-(i/,7) i:{n)~>{N) se{n) 

We thus have 

(LHS of dMl)) 

= E E E E 

i:(n)-i.(Af> k:'X-i.(A'> j:'I-i-(n) l/>:'I-*(n) 
s.t.jop— j onc-to-onc 

Il{7r ^ v o '0}]l{j =70 T/ijIlIk = i o T/jjlIk op = io9oip} 

n 2;^M,7M(i(s)>i(^(s))) 

se{n)\^{i) 

= E E n iHs)Ms)iKs)M0{s))) 

Jp:i^{n) i■.{n)-^{N) se(n)\i/'('I) 

ono-to-one s.t. ioi/iop^ioeoi/i 

s.t. UOIp—TT 

and 'yoxjjop—'yo'tp 

= Yi Yi n ^^(^),7(s)(K5), 1(^(5))) 

rGSplice^(:/,7) i:(n)-^(iV) se(n) 

s.t. i(0r(s))— i(s) s.t. r(s) — s 

for s £ (n) s.t. t(s) ^ s 

= (RHS of dMl), 
which finishes the proof. D 

4.2. Concluding arguments. 

4.2.1. Rewrite of the right side of ©. Let E^p be the A^-by-A^ GUE matrix to 

which the matrix X^- evaluates when all the variables Xij{a,l3) are evaluated at 
i.i.d. standard normal random variables. Then by Theorem 12.4.51 combined with 
Proposition I4.1.4[ with J2% denoting summation over trees T spanning (k) and 
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picking the random matrices wtj independent of all the GUE matrices S^ , we 
have the formula 

, , / expression under the limit on the right 

^ ' \ side of © muhiplied by iV2+"/2-fe 



Ee 



tr Polye 



>7i'^i-i- 



X,=Z.r(<^(i),7(i)) 
for i £ (n) , 



for {i.j) £ (k) X (n) 



Etr Arbe, 



7,1>' 



for (i. j) S (fc) X (n) 



4.2.2. Passage to the limit. Divide through ([27]) by iv2+"/2"'''. Consider the limit 
as A^ — > oo on extreme left and extreme right. On the one hand, by Theorem 1 1.1 1 on 
the left, we recover the quantity |Mapg(0,7)|. On the other hand, by Voiculescu's 
result [32, Thm. 2.2] on the right, we recover the right side of p7[) . The proof of 
Theorem 12.3.51 is now complete. D 

Remark 4.2.3. The preceding argument provides a proof of the existence of the 
limit on the right side of ^ independent of the usual diagram-based approach. 

Remark 4.2.4. Since relation (|27|) already holds exactly for each fixed N, the meth- 
ods introduced in this paper seem to be susceptible to higher genus generalization. 
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